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Abstract 



To have an analytic understanding of the higher-twist structure functions, 
^ I we calculate twist four longitudinal structure function for a positronium-like 

bound state in weak coupling light-front QED. We find that in the weakly 
coupled system, the fermionic part of Fl is related to the kinetic energy of 
the fermions and not to the interaction. We verify a previously proposed 
sum rule in this limit, which in this case reduces to a relation connecting the 



if) , kinetic and the potential energies to the binding energy of positronium. Using 

I the analytic form of the wave function of positronium in this limit, we show 

■<!:;;j- ' that the constituent counting rule does not hold for x — > 1. The twist four 

■ Fl in this limit is similar in form to a widely used phenomenological ansatz. 

fH , Keywords: Light-front Hamiltonian, hound states, twist four, structure function 

p • 1. Introduction 

^ ■ 

Higher twist or power suppressed contributions to deep inelastic scattering structure 
functions involve non-trivial non-perturbative information about the structure of hadrons. 
^ I Recent experimental results indicate that these higher twist effects play an important role in 
d ' the kinematical range of the SLAC experiments and it is important to have a clear physical 
picture of them. Light-front Hamiltonian QCD offers a theoretical tool to investigate the 
deep inelastic scattering structure functions. This is based on physical intuitions and at 
the same time employs well defined field theoretical calculational techniques. The structure 
functions are expressed as the Fourier transform of the matrix elements of light-front bilocal 
currents. Fock space expansion of the target state allows us to express these in terms of 
light-front multiparton wave functions. The interesting aspect of this formulation is that 
both perturbative and non-perturbative issues can be addressed within the same framework 
IjlJ. Recently, the twist four part of the longitudinal structure function Fl and the transverse 
polarized structure function qt, which is a twist three contribution have been analyzed in 
this approach and various interesting issues associated with them have been addressed 
The structure functions can be calculated once the light-front bound state wave functions 
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are known. However, the actual calculation for a QCD bound state in 3 + 1 dimension is 
highly complicated and it requires the recently developed similarity renormalization tech- 
niques [0 for the light-front QCD Hamiltonian. The spectra of the heavy quark bound 
states like charmonium and bottomonium have been investigated using this technique [^. 
The analytic form of the wave function is not obtained so far. This is because, the similar- 
ity renormalization technique generates an additional confining interaction in the effective 
Hamiltonian which makes it impossible to solve the effective Hamiltonian analytically even 
in the leading order in bound state Hamiltonian perturbation theory. In this work, we have 
performed a much simpler but quite interesting analysis. We have calculated the twist-four 
part of the longitudinal structure function for a positronium-like bound state in light-front 
QED in the weak coupling limit. In fact, in this limit, QCD results are not expected to 
differ too much from the QED results. As a result, this analysis is important since it tests 
and illustrates the approach in QCD. The advantage here is that, in the leading order in 
bound state perturbation theory, the bound state equation can be solved analytically and 
the wave function is known. This allows an analytic understanding of the problem. 

In the weak coupling limit, we show that our previously proposed sum rule reduces to 
a relation connecting the kinetic and potential energies to the binding energy of positronium. 
Using the analytic wave function of positronium in this limit, we show that the structure 
functions fall faster than that predicted by constituent counting rule near x —>■ 1. Also, we 
find that the twist four in the weak coupling limit has a form which is similar to a widely 
used phenomenological formula, which connects the twist four distribution to the twist two 
distribution. 

2. Twist Four Longitudinal Structure Function 

We consider a positronium-like bound state | P) given by 
dk-^ d k-^ r c/A^2 d k'2 



MP \ki,(r,; k2, a2)p{{2T,YP+5\P - k, - k^Pik^, ai)d\k2, as) | 0) 
^ — ^ r dk^ d k-^ r dk^ d k^ f dk^ d k^ 
a^7^,xj ^2{2Tifkt ' ^2{2Tifkt ' ^2{2^fkt 

MP I fci, ai; A;2, a^, k^, Xs)^2{27r)^P+6''{P - k^ - k^ - k^) 

fet(fci,ai)dt(A;2,a2)at(A;3,A3)|0). (1) 

Here 02 is the probability amplitude to find an electron and positron in the positronium, 
03 is the probability amplitude to find an electron, positron and a photon in the positronium. 
We consider upto three particle sector. We introduce Jacobi momenta (x, k-^) and the boost 
invariant amplitudes, -02 and V'a HI- 

We calculate Fl^^ in light-front QED in light front gauge, = using weak coupling 
approximation. This includes the effect of dynamical photon. It can be shown that, for a 
weak coupling theory the results are the same as obtained using a non-relativistic approx- 
imation. However, the entire calculation is fully relativistic and exact in the leading order 
in bound state perturbation theory 0. 
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The twist-4 part of the fermionic component of the longitudinal structure function is 
^iven by, 



(2) 



Ml 



1 x2fP^ 



Q2 2% 



dy- (T^P^y-- {P I ^(y-)7-^(0) - ^(0)7-^(1/-) I P), 



(3) 



Mo 



{P^f 1 x2(P+) 



(P+)2 g2 27r 



£ / rfy- e-5^^J^-- (P I V^(y-)7+^(0) - V^(0)7+V^(r) I P)- (4) 



We shall take the mass of the state | P) to be M and the electron and positron mass 
to be m. In the weak coupling (non-relativistic) limit the helicity dependence of the wave 
function factorizes away, so it is sufficient to consider one helicity sector. Here we shall take 
the two particle state with helicities ai and o"2 up. 

For a positronium state | P) given by Eq. (|T]) we obtain. 



pr=4 , 



{M,)diag + (M2) diag 

+ ^ ^ y dxid'^K,^ J dx2d'^K2 \ 1^3 f (('^i')^ + rn'^)S{x — xi) + 
{{k^Y + m'^)5{x - X2) . 
The off-diagonal contributions to P£y-j(x) comes from Mi alone. 

4 e2 



(Mi)„ 



where 



Ml" 



ff-diag 



1 



(xi - y)'^E 



2(27r)3 



rfxirf^K^ / dyd^K^ Mf + Ml' + + Mf 



Xi 



8{x - xi) - - ^)] I ^s'^^'l^i- 4) 1', 



(5) 



(6) 



(7) 



M° 



i?(x — y)2 



- (/t]^)^5(x - xi) + {K-^f5{x - y) 
(V^r'^^(xi,/s:^)V^-^-^(y,«:^) + /i.c.), 



(8) 



Ml 



1 r 2(4)^y 
(2/ - xi)2E' Kl - Xi) 



(5(1 — X — xi) 



2(k:^)2(1-xi) 



6{x-l+y) 



(xi, /ti 



(9) 
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{K^f5{x + y-l)- {K^f5{l - X - xi) 



(10) 



In this calculation we have taken all operators to be normal ordered. Also, in the Eq. (|^- 
(p!0|) we have neglected all mass terms in the vertex, since these terms are suppressed in the 
non-relativistic limit. The energy denominators are given by. 





(4- 




y I - xi 


Xi 


-y 








xi y 


iy- 


Xi) 



E = M' 
E' = 

We define the twist four longitudinal photon structure function as 

1 xP^ 



(11) 



Q2 2n 



dy- e-^^+^ 

(p+)2 



1 



{P I {-)F^\y-)F^{0) + -g^-F''^{y-)F,M \ P) 



{P\F+\y-)F^{0)\P) + {y--0) 



(12) 



Fl^Jj{x) has both diagonal and off-diagonal parts. We take all operators to be normal 
ordered and we get. 
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X diag 

4 4e2 



dxid Ki / dyd k \ ip^ 



(1 - xi - y) 



■5(1 — xi — x — y) 



dxid'^K^ J dyd^K'^ip*2{xi. 



{xi -y) 



-6{xi -x-y). 



(13) 



The second term in the right hand side is the contribution of the instantaneous interaction. 
The off-diagonal contribution is. 



nT=4/ 



^9) off -diag 



G1 + G2 



(14) 



where. 



^1= 



4 e2 



■ j dxid^Ki j dyd^K^ 



X 



(-4 



g2 2(2^)3 7 W " E{xi-yy 

(4 - K^f 2(4)2 2{K^f 



{xi - y) 



+ 



Xi 



2(4-^:^)2^ (4)2 



±\2 



Xi-y) (1-Xi) (1-?/)' 



*0"lO"2 



(15) 



(- 



g2 2(27r)3 J ' ^ J ^ E' 
4(4 - K^f 2{K^f 2(4)^ 



2(4 - 



(^r'^'la^i, 't^)^2'"H2/, '«^) + h.c.)]6ixi + x-y). 



Where E and E' are given by Eq. (|rTj). 
From these expressions, we calculate, 



1 F£g(x) + Fl^f^jx) 

X 



2 



dx = — f dxd^K^4j2*H^^ + 7 
J X 1 



X 



+ 



4 

g^2(2^ 

4 4e2 



J dxd^K-^ j dyd^q^[Ai + A2 + Bi + B2] 
J dxd^K'^ J dyd^q^- 



2(27r)3 

i)l{x,ti^)i)2{y,q 



where, 



1 1 



E{x-y) 



A, 



1 1 



E' {y - x) 



1 1 

E{x-y) 



Bo 



1 1 



E' {y — X 

where Vi,V2,V'i and V'2 are given by. 



y'2^2"^'^^(x, K^)V'r'"'(?/, + ^'il V^2^'''(2^, 



(x — {x — y)^ x^ 



(k^)2(1-2x) ^ iq^)\2y-l) 2{k^ - q^f 



x{l-x){x-y) y{l-y){x-y) (x - yf 



(24) 



In these expressions we have kept only those terms in the vertex which survive in the 
non-relativistic hmit. The hehcities ai and a"2 are both up and we have neglected all mass 
terms in the vertex since in this limit they are suppressed. Also in the weak coupling (non- 
relativistic) limit, we consider only photon exchange interactions and the off-diagonal terms 
proportional to | 7/^2 I originating from self energy effects can be neglected. 

Now, from the expression of -03 (see 0) we can write, 

„2 



I dxd\^ J dyd^q^l^B, + -is2 



where Bi and B2 are given earlier. Using this, one can write the second term in the right 
hand side of Eq. (p^Tf) as. 



g2 2(27r)^ 



-J dxd^K^I dyd\^[^B, + ^B 



+ 



+ m — K — q 



X y il-x-y) 

Considering the fact that the total energy is conserved, one can write this as. 



g2 2(27r) 



j dxd^K^ j dyd^q^{Bi + B2) 



(26) 



(27) 



So we get. 



1 Fl^^{x) + F, 



T=4, 



X] 



-dx- 



X 



dxd^ K'^ il)2* ^2 



±\2 



+ 



X 



X 



j dxd^n^ j dyd^q-^[Ai + A2 
dxd^K'^ I dyd^q 



4 Ae 



2(27r) 



{x - yf 



(28) 



Also, if we denote 



x{l—x) 



by M^, then in the non-relativistic limit, it can be shown that, 



— Mq ~ O(e^). So we neglect this difference in the energy denominators and replace the 
bound state mass in E and E' by = The energy denominators then become, 

E 



x(l—x) 



X 



x-y 



1 



,m 



{x — y) X 



[{-y{x-yY + {K' 



(29) 



E' 



+ m 



[K-^ - q^ 



1 — X 
1 



ix-y)h 



m 



X 



1-y x-y 
^\x - yf + {n^ - q^f\ 



(30) 
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We get, in this limit, 



i 



1 Fl='{x) 



dx 



X 



4 r 

^ J dxd^K I ^2 



Q 

4 2e2 



xfl — x) 



■ j dyd^q^%l)*2{x, k^,1-x, -K^)'4)2{y, q^,l-y, -q^ 



1 



(31) 



Here in the weak couphng hmit we have omitted the spin indices. 
The Fermionic part of the Hamiltonian density is given by, 

0+- = i^j-d+i; = 2^+^^ \a^.{id^ + gA^) + ^^m] la^.Ud^ + gA^) + j^mW 
J I J I J 

The gauge bosonic part of the Hamiltonian density is given by 



(32) 



+ (^) 2(^+)V+ (^) 2{r)'r. 



d+ 



d+ 



(33) 



For a 



The fermionic part of the longitudinal momentum density is given by, 9j 
and the gauge bosonic part of the longitudinal momentum density, 6'^^ 
positronium-like bound state, we calculate the matrix element of 9j~ and 9'^~ . The matrix 
elements have both diagonal and off diagonal contribution. The diagonal contribution to 
the matrix element from the fermionic and the gauge bosonic part is given by. 



(p±)2 

{p 1 0+- I p) - )jr-^{p I e++ I P) 



(p^y 



2 j dxid^K^'il;*ijj2{ 



+ 



2 J dxid^K^ j dX2d'^K2lp^tl>3{ 



diag •' Xi [L 

{>^ir , (4) , i-4 



Xi 



■} + 



+ 



8e' 



2(27r): 



Xi X2 (1 - Xi - X2) 

j dxid^Ki j dyd'^K'^ip2{xi^ Ki 

^2(2/, K^) 7 ^-T^ 

{x - yy 



} 



(34) 



where 9+- = 9j- + 9+' 

The off-diagonal part can be written as. 



(p±)2 

{P I 9+- I P) - f^{P I 9++ I P) 



V1 + V2 



(35) 



off-diag 
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where, 



Vi 



2e^ 



2(27r)= 



dxd^ 



K 



dyd^q 



2Vi\ i/j^'''%x,K^) 



1 1 

E{x-y) 



(36) 



2e^ 



dxd^K'^ / dyd^q 



1 



2(27r)3 7 J " E'{y-x)^ ' ^ ' ^ 

+y'2(V'r''Ha;, K^,l-x, -K^)ilj^^'^'{y, q^) + h.c.) 



(37) 



where the expressions for Vi,V2,V'i,V'2, E and E' are given earher. As before, we have 
taken the two particle state with both ai ,a2 up. 

Considering only the photon exchange interactions and putting = Mq in the energy 
denominators as before, one obtains in the non-relativistic limit for a weak coupling theory. 



{p I e+- 1 p) 



„ 2e2 



(P I I P) =2 dxd'K^ \ i); 



)2 + 



x{l — x) 



2(27r)= 



g ij*2{x,K ,i-x,-K )ip2{y,q ,i-y,-q ) 



X 



1 



+ 



m 



+ 



"m.\2i 



X - yY 



:i-x) 



+ 



'{x-yf 



(38) 



We can see that the right hand side of the above equation, which is nothing but Coulomb 
interaction, has exactly the same form as the interaction part of /q ^ ^ . So, the interac- 
tion part of /q ^ ^ dx can be related to the expectation value of the Coulomb interaction. 



Since k and q are small in the non-relativistic limit, all 



and {q^y dependence in 



the numerator of the interaction terms are neglected compared to the m dependent terms. 



However, the term proportional to 



cannot be neglected because both x and y are 



almost equal and this term cancels the instantaneous interaction in the non-relativistic limit. 
Both of these terms originate from PL(g) and one can see that only the gauge bosonic part 
of the longitudinal structure function is important for the Coulomb interaction in the weak 
coupling limit. The terms in the energy denominators combine with the other terms to 
give non-vanishing contribution. The fermionic part only gives contribution to the kinetic 
energy of the fermions. Reminding oneself that we are working in the light-front gauge and 
not in the Coulomb gauge, this is a manifestation of the gauge invariance of the separation 
of the Hamiltonian density into a fermionic and a gauge bosonic part. 

We introduce the three vector p as, p = {k, k^) where is defined through a coordinate 
transformation from x G [0, 1] to G [— oo, oo] by, x = | H ; , .'^^ , =. We introduce 



the bound state wave function (^{p) which is normalized as, / d^p(p*{p)(p{p) = 1. The bound 
state equation for the positronium in the weak coupling limit can be written as [§, 



-4((p)2 + m2) 0(p) 



2e' 



2(2 



d^p4>{p) 



Am 



{p-p') 



/\2 



(39) 
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We write Eq. (gH]) as, 



{P I r- I P) - |^(P I I P)]=2j d'p\ m \h[{p)' + m'] 



We can now see the more familiar form of the Coulomb interaction. Multiplying the bound 
state equation Eq. (^) by (/)*{p) and integrating we get, 

M' = l £p\ m Mp)' + m'] - I d'pj d'p'^*(p)^(p')-0^. (41) 

Hence from Eq. (|3lD , (|38| ) and ( PD it can be seen that the sum rule [|| is satisfied in the 
weak coupling limit for a positronium target in light-front QED and it can be written as, 

1 Fr=^(x) 2 r (P-^V 1 

^ A [(P I r- I P> - I I P}] ^ 4^. (42) 

In the non-relativistic limit for a weak coupling theory, = 4m^ + AmBf. where B^. is the 
binding energy of positronium. 
From Eq. (|4l|) we obtain, 

= ^(P) - ^ /''V/ (43) 

The first term in the right hand side is the kinetic energy with y being the reduced mass 
of the two body system and the second term is the expectation value of the Coulomb 
interaction. So we see that in the weak coupling limit, the sum rule reduces to a relation 
connecting the kinetic and potential energies to the binding energy. 

3. Fl^^ for the ground state of positronium 

The bound state equation Eq. (^) can be analytically solved for QED, which is the 
primary motivation for studying QED. The ground state wave function of positronium is 
given by, 

0!^,Se,v(^' ^' = (l^y{P)^S,,S^S^,,s' (44) 

where Se and Se' label the spin quantum numbers of the electron and positron respectively 
and V denotes all the other quantum numbers, p = n,l,m correspond with the standard 
non-relativistic quantum numbers of hydrogen. The spin part factorizes out and the wave 
function is normalized to 1. The wave function is given by, 

where e„ = ^ and Y,^{Qp) = y„;„j(f2) are Hyperspherical harmonics. Here <| m |< 
Kn-1. 
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For Is state of positronium, we have, 1^1,0,0 
wave function can be written as, 



-X=. In terms of x and k"*-, the Is state 



4(ei) 



ei)2 — + 



1 (K^)2+m2 
4 x{l-x) 



l2 



(46) 



which agrees with [^] for non-relativistic — |- 

The weak couphng hmit contributions to the structure functions F2{x) and Fl^^{x) can 
be directly evaluated using this wave function. 

|2 



AxHl 



XI 



K 



(47) 



where A = 4.78 x 10'^^{Mevf. 
introducing a cutoff A. We get, 



m'^[{l - 2xy + a^x{l - x)] 
The integral is convergent and can be evaluated analytically 



F-iix) = 28.15 X 10 



-11 



X 



(1 -2x)2 + a2x(l 



X] 



3 • 



(48) 



F2{x) is very sharply peaked at a; = |. 

The twist four longitudinal structure function is given by, 

\2 



X 



4 

4 



X 1 



x) 



)vri 



x) 



m2[(l - 2.t)2 + a'^x{l - x)] + {K^f 



4 • 



(49) 



Here we have considered only the (k-'-)^ dependent part, the integral of which is directly 
connected to the kinetic energy of the electron-positron pair. The above form is similar to 
the widely used phenomenological ansatz for the twist four distribution [|T^ but not exactly 
the same. Evaluating the integral analytically, we get. 



Fr(x) 



X 



14.7 X 10 



-11 



1 



X (1 — x) 



(1 - 2xy + a^x{l-x) 



2 • 



(50) 
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is maximum at x = ^ and it is also sharply peaked. The structure function calculated 
above (|48|) falls faster than that expected from constituent counting rule near x 
This is because the non-relativistic wave function behaves differently in the asymptotic 
region. It can be seen that it behaves as 



for large k whereas the relativistic wave 
function behaves as vprw [0- The asymptotic behavior of the form factor of positronium 



F(Q ) can be calculated analytically using the above wave function. It can be shown that 



F{Q^) ~ -Ai when is large, as expected for positronium in the non-relativistic limit |]T2 
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This means that it docs not have a monopole behavior in this hmit and one cannot expect 
the constituent counting rule to hold. 

To summarize, in this work, we have investigated the twist-four longitudinal structure 
function for a positronium-like bound state in light-front QED in the weak coupling limit. 
In this limit, we get expressions that look similar to the familiar non-relativistic expressions, 
but the entire calculation is fully relativistic in the leading order in light-front bound state 
perturbation theory. We have explicitly verified a sum rule for -F^^^ that we previously 
proposed. It is worth mentioning here that in the non-perturbative context, we have investi- 
gated before the twist four longitudinal structure function for a meson in 1-1-1 dimension and 
showed that the sum rule is satisfied using t'Hooft's equation. In this work, we have shown 
that in the weak coupling limit, the sum rule reduces to a relation connecting the kinetic and 
the potential energies to the binding energy of positronium. We have also shown that, in 
this limit, the fermionic part of FJ^""^ contributes only to the kinetic energy of the fermions 
and not to the interactions. The twist four in this hmit looks similar to a much used 
phenomenological ansatz, however, here we get the result directly from field theory. The 
structure function falls faster than that predicted by constituent counting rule, as expected. 

The twist four part of the longitudinal structure function is important since it is the 
leading non-perturbative contribution to F^. The leading twist contribution to Fj, is per- 
turbative, in contrast to the case of F2. This analysis for a bound state in weak-coupling 
light-front QED is quite interesting since it gives an idea of what goes in such a calculation 
in light-front QCD. Similarity renormalization upto O(e^) does not produce any additional 
interaction in the efi^ective QED Hamiltonian and in the weak coupling limit, one can work 
with the canonical Hamiltonian. The overall computational framework in QCD is the same 
and this analysis in light-front QED allows an analytic understanding of the problem. 

I would hke to thank Prof. A. Harindranath for suggesting this problem and also for 
many illuminating discussions. 
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